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Outline 

Part  I  -  Snow  mechanics 

-  Background 

-  Experimental  procedure 

•  Tribometer  for  indentation,  piowing,  siiding  tests 

•  3D  X-Ray  Microtomography  for  microstructure 

-  Numerical  modeling  procedure 

-  Typical  results  (indentation,  plowing,  compression 
tension,  penetration) 

Part  II  -  Vehicle-snow  interaction 

-  Alaska  Instrumented  Vehicle  and  profilometer 

-  Validation  of  models 
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Background: 

Characteristics  of  (Geometric)  Snow  Models 


•  Multi-scale  in  nature: 

-  um  scale  at  the  sub-grain  level  (microscale) 

-  mm  scale  at  the  grain  level  (mesoscale) 

-  cm  scale  at  the  terrain  level  (macroscale) 

•  Stochastic  in  nature: 

-  Stochastic  models  at  each  scale  (e.g.,  Gaussian  Random 
Field  at  the  mesoscale,  semi-variogram  at  the  macroscale) 

-  Key  challenge: 

•  Integrate  (‘patch’)  models  at  different  scales 
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Background: 

Ipdentation,  plowing  and  sliding 
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(c)  top  view 


Resultant  Forces  due  to  Sinkage/Ploughing  and 
Longitudinal/Lateral  Slips 

Motion  Resistance,  Shear  Force  and  Drawbar 
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Background: 

Needs 

Microstructure  (uncertainty)  effect  not  assessed 

Need  better  understanding  of  deformation  and  faiiure 
mechanisms 

Littie  work  done  in  piowing  and  siiding 
Size  effect  not  understood 
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Background: 

Goals  and  Approaches 

Goals: 

-  Develop  models  for  the  mechanical  properties  of  different 
types  of  snow 

-  Quantify  the  associated  uncertainties  and  understand  the 
sources  of  uncertainties 

Approaches: 

-  Experimental: 

•  Microscale  tests  using  microtribometer 

•  Microstructural  statistics  using  microCT  scanner 

-  Numerical: 

•  Microscale  simulations  using  a  meshless  method  with  appropriate 
constitutive  laws 

-  Semi-analytical: 

•  Continuum  mechanics  based  stochastic  models  incorporating 
microstructural  information 
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Experimental  Procedure 


Collection  and  storage  of  snow 

-  February  to  March,  2009,  Tanana  River,  Fairbanks,  Alaska 

-  Fine-grained  just  underneath  the  surface 

-  Coarse-grained  about  20  cm  from  surface 

-  Snow  temperature  ~-6  C 

-  Stored  in  a  freezer  ~-25  C 

Microtribometer- 

-  Temperature  ~-10C 

-  Pin  sizes  (1/8”,  1/4”,  3/8”,  1/2”) 

-  Force  or  velocity  control 

-  Multiple  steps  and  modes  (indentation,  pin-on-disk  etc.) 
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Experimental  Procedure: 
tribometer  setup 


Environment 
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Pin-on-disc  setup 

Approved  for  public  release 


Experimental  Procedure: 
Skyscan  1172  Microtomography 
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Experimental  Procedure: 
Snow  Sample  Holder 
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Experimental  Procedure: 
Grey-level  Cross-Sectional  Image 
Sieved  Snow  <  1  mm  Grain  Size 


Brighter  pixels 
represent  ice 


7.344  mm  by  7.344  mm 
Resolution:  1225  by  1225,  Pixel  size:  6  micron 
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Experimental  Procedure: 
Grey- Level  Histogram 

Histogram 
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Experimental  Procedure: 
Segmentation 


Black  is  ice 


grey-level 


binarized  image 


UNCLAS:  Dist  A.  Approved  for  public  release 


Experimental  Procedure: 
Removal  of  Unconnected  Parts 


speckles 


Binarized  image 


Remove  speckles 
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Experimental  Procedure: 

3-D  Visualization  of  a  Cube  of  Snow  Microstructure 

Side  Length  =  3.618  mm 
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Experimental  Procedure: 

Extract  Statistical  Information  from  Images 


Porosity  & 
(pore  volume 
fraction) 


Probability  that  two 
points  a  distance  r 
apart  will  lie  in  pore 
space 


Two-point  probability 
function 
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Numerical  Modeling: 

Generalized  Interpolation  Materiai  Point  (GIMP)  method 

(1/2) 

•  Geometry  from  CT  images 

-  148x148x148  voxels  (48  um  resolution); 

7.1mmx7.1mmx7.1mm 

-  Each  voxel  (ice)  is  mapped  to  a  material  point  (particle) 

-  ~1  million  particles 

•  Boundary  conditions 

-  Periodic  on  the  sides  (for  indentation) 

-  Frictionless 

-  Speed  of  indentation  is  71  mm/sec 

•  Indenters 

-  1/16”,  1/8”,  1/4” 
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Generalized  Interpolation  Material  Point  (GIMP)  Method 

(2/2) 

•  Software:  parallel  code  Uintah  installed  on  a  Sun 
cluster  at  Arctic  Region  Supercomputing  Center 

•  Constitutive  law  used  for  ice  particles 

-  Elastic-brittle  [cf.  Johnson  &  Schneebeli  (1999),  Marshall 
and  Johnson  (2009)] 

•  Failure  according  to  maximum  tensile  stress 

•  Post  failure 

-  Stress  set  to  zero  if  mean  stress  is  tensile 

-  Stress  set  to  mean  stress  if  compressive 

•  Algorithm 

-  Dynamic,  explicit 
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Tests  and  Simulations 


Tests 

-  Compression 

-  Indentation 

-  Plowing 

-  Sliding  on  compacted  snow  (future  work) 

-  Penetration  (future  work) 

Simulations 

-  Compression  and  Tension 

-  Indentation 

-  Plowing 

-  Sliding  (future  work) 

-  Penetration 

-  Triaxial  tests 
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Typical  Results:  Indentation  tests  for  fine  snow 


Fine-grained  snow  depth=18mm,  speed=5mm/sec 
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Microstructure  after  Indentation  Tests  via 

MicroCT 


Top  View 


Fine-grained  snow: 

initiai  density:  ~290  kg/m'^3 
Finai  density:~590  kg/m'^3 


Side  View 
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Typical  Indentation  Simulation  Results 


Fine-grained  snow,  71  mm/sec  1 /4-inch  indenter 
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Normalized  cumulative  damage 


Normalized  cumulative  damage 


Typical  Indentation  Simulation  Results 

Cumulative  damage 


Fine-grained  snow,  71  mm/sec 
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Failed  Particles  from  Indentation 

Simulation 


DB:  index.xml 


Morijun  15  1 1:44;07  2009 
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Characteristics  of  Indentation  Test  Curves 


Fine-grained  snow,  5mm/sec 
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Background: 

Indentation  modeling  using  continuum 

mechanics 


Pressure 


Three  zones: 

I:  Elastic 

II:  Hardening  (via  cavity 
expansion  theory  and 
Drucker-Prager  criterion) 

III:  Densification  (via 
upper  bound  theory  and 
Drucker-Prager  criterion) 


*J.H.  Lee,  J.  of  Terramechanics  (2009) 
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Potential  Deformation  Mechanisms 


A:  Upper  ‘yield’  point 
(inelastic  due  to 
damage) 

B:  Lower  ‘yield’  point 
-  OAB:  Initial  yield  zone 

B-C:  Hardening 
(additional  damage) 

C:  Plateau  stress 


C-D:  Compaction  (little 
additional  damage) 

D-E:  Densification 
(pressure  bulb  hits 
bottom) 

UNCLAS:  Dist  A.  Approved  for  public  release 
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Results:  Plowing  tests 


Fine-grained,  1 /4-inch 
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Load  (N) 


Results:  Snow  Penetration  Simulations 
(45  deg  inclusion  angle)* 


*Lee  et  al.,  Proceedings  of  ISTVS  2009 
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Results:  Typical  Penetration  Geometry 
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Results:  Strengths  from  Inversion  of  Penetration 

Signals 
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Part  II: 

Vehicle-Snow  Interaction 


An  instrumented  vehicle  (Alaska 
Instrumented  Vehicle)  to  collect  data  about 
vehicle  and  wheel  states 

A  vehicle-mounted  profilometer  to  measure 
terrain  topology 

Equipment  to  obtain  microstructure  and 
mechanical  properties  of  snow 
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Alaska  Instrumented  Vehicle 


2008  Jeep  Commander  (with  ESP) 

Vehicle  states: 

-  Longitudinal  slip  (via  wheel  longitudinal  speed  and 
wheel  angular  speed  from  ESP) 

-  Vehicle  speed,  sideslip,  wheel  slip  angle,  yaw, 
pitch  and  roll  (VBOX  II  SX  ?+  ESP) 

-  Wheel  forces  and  moments 

•  Kistler’s  wheel-force  transducers  (a  set  of  4) 

Validation  on  pavement  first 
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Terrain  Profiling 

Vehicle-mounted  profilometer  (Kern  and 
Ferris,  2007) 

-  Inertial  navigation  system  (INS)  to  determine  the 
position  and  orientation  of  the  vehicle 

•  Differential  GPS  system 

•  Inertial  measurement  unit  (IMU)  -  gyros  and 
accelerometers  for  orientation  and  position 

-  Scanning  laser  for  profiling 

-  4-meter  wide  scan  (claimed  accuracy  of  vertical 
measurements  0.7-1 .0  mm) 

-  Claimed  horizontal  precision  is  1mm  for  short- 
distance  traveled 
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Measurements  Needed 

Depth  of  snow  cover  ~5  cm  -  30  cm 

Snow  density  and  in-situ  compressive 
strength 

Mechanical  properties  and  microstructure  by 
collecting  and  transporting  select  samples 
from  field  to  lab 

Vehicle  and  wheel  states 
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Tentative  Test  Protocols: 

Before  Vehicle  Travel 

Select  areas  for  types  of  snow  -  (dry,  wet, 
windblown  etc.),  depth  of  snow,  strength  of 
snow  -  with  enough  room  to  maneuver  the 
two  vehicles  (AIV  and  profilometer) 

Measure  snow  depth  by  profiling  ground 
twice  -  with  and  without  snow  (winter  first, 
summer  later) 

Measure  snow  properties  along  the  intended 
path  before  vehicle  travel 
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Tentative  Test  Protocols 


•  Passes: 

-  Single  pass:  rut  created  by  front  wheels  not 
traveled  by  rear  wheels  for  virgin  snow 

-  Multiple  passes  for  compacted  snow 

•  After  vehicle  travel: 

-  Measure  sinkage  (3D)  using  profilometer 

-  Measure  deformed  mechanical  properties  of  snow 

•  Maneuvers: 

-  Combination  of  driven  and  driving  wheels 

-  Longitudinal  and  lateral  motions 

-  Effects  of  ESP 
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Development  and  Validation  of  Models  for 

Virtual  Proving  Ground 

Development  of  stochastic  terrain  models 

Improvement  of  indentation  model  (J.  Lee, 
2009) 

Validation  of  stochastic  tire-snow  interaction 
model  for  combined  slip  (Li  et  al.,  2009) 

Validation  of  finite  element  tire-snow 
interaction  model  for  combined  slip  (J.  Lee, 
under  review) 

Validation  of  time-dependent  tire-snow 
interaction  model  for  combined  slip  (Lee  and 
Liu,  2006) 
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